Matter wave interference using two-level atoms and resonant optical fields 
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A theory of matter wave interference is developed in which resonant, standing wave optical fields 
interact with an ensemble of two- level atoms. If effects related to the recoil the atoms undergo 
on absorbing or emitting radiation are neglected, the total atomic density is spatially uniform. 
However, when recoil effects are included, spatial modulation of the atomic density can occur for 
times that are greater than or comparable with the inverse recoil frequency, w" 1 . In this regime, 
, the atoms exhibit matter-wave interference that can be used as the basis of a matter wave atom 

■ interferometer. Two specific atom field geometries are considered. In the first, atoms characterized 

by a homogeneous velocity distribution are subjected to a single radiation pulse. The pulse excites 
the atoms which then decay back to the lower state. The spatial modulation of the total atomic 
density is calculated as a function of t, where t is the time following the pulse. In contrast to 
the normal Talbot effect, the spatially modulated density is not a periodic function of t, owing 
to spontaneous emission; however, after a sufficiently long time, the contribution from spontaneous 
O |' processes no longer plays a role and the Talbot periodicity is restored. In the second atom-field 

geometry, there are two pulses separated by an interval T. The atomic velocity distribution in 
this case is assumed to be inhomogeneously broadened. Owing to the inhomogeneous broadening, 
one finds a nonvanishing spatial modulation of the density only at specific "echo times" following 
the second pulse. In contrast to the normal Talbot-Lau effect, the spatially modulated density is 
not a periodic function of T, owing to spontaneous emission; however, for sufficiently long time, the 
contribution from spontaneous processes no longer plays a role and the Talbot periodicity is restored. 
• i-H . The structure of the spatially modulated density in the vicinity of the echo times is studied, and is 

found to mirror the atomic density following the first pulse. With a suitable choice of observation 
time and field strengths, the spatially modulated atomic density serves as an indirect probe of the 
Q ■ distribution of spontaneously emitted radiation. 
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Atom interferometry has emerged as an important new research area in the past fifteen years Q . Thermal atoms 
offer unique properties as the working element of an interferometer, owing to their small de Broglie wavelength. 
Atom interferometry has already led to impressive results in precision measurements of h atomic polarizability 
[3|, rotation rates [|]J|] and the acceleration of gravity ||. Applications to nanolithography have also been proposed 
[7|. Despite the advances in atom interferometry, there remains some controversy as to what constitutes an atom 
interferometer. We have discussed previously j|] the distinction between classical and quantum scattering for atom 
interferometers in which microfabricated gratings are used as the "optical" elements in the interferometer. In this 
paper, we extend this discussion to atom interferometers that employ optical fields to modify the atomic state wave 
function. 

In our effort to classify atom interferometers, it will prove useful to summarize the conclusions we reached in 
analyzing the scattering of atoms by microfabricated gratings ||. As a model system, we consider the scattering of 
an atomic beam by two microfabricated gratings separated by a distance L. The internal structure of the atoms is 
irrelevant here, since the atoms are in their ground state and the scattering by the gratings is assumed to be the same 
for any and all ground state sublevels - there is no internal state coherence. For an angular divergence of the atomic 
beam, 9b, satisfying 9b S> d/L, where d is the grating period, the atomic density is uniform in the region between the 
gratings, but can exhibit "fringes" beyond the second grating. In particular, at a distance 2L from the first grating, 
fringes having period d are produced. As such, this device has the appearance of an interferometer. Moreover, it can 
be used to measure inertial effects |J , since the fringes are displaced if the atoms accelerate between the interaction 
zones. Further evidence that this device constitutes an interferometer is provided by a theoretical description of the 
scattering in which the microfabricated gratings create coherence between transverse momentum components of the 
incoming atomic beam, taken as a plane wave. Despite all these features, it is incorrect to conclude that such a system 
legitimately qualifies as an atom interferometer. 
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To clarify matters, it is best to back up a bit and consider the scattering of an atomic beam by a single grating. At 
a distance L from the grating, the nature of the scattering is determined by the so called Talbot length |1CJ], defined 

as 

Lt = 2d 2 /XdB, 

where XdB — h/muo is the atomic de Broglie wavelength (m is the atomic mass and uq is the longitudinal atomic 
speed). For L <C Lt, diffraction effects are unimportant and one is in the classical scattering limit. For a highly 
collimated beam (9b <C d/L), the pattern produced on a screen at a distance L <C Lt from the microfabricated grating 
is simply the shadow of the grating. This is clearly a classical effect. 

Of course, the classical scattering limit can be obtained from a calculation in which the center-of-mass motion 
is quantized. If such a fully quantum calculation is carried out, one finds that, immediately following the grating, 
the atomic state wave function contains coherence between different transverse momentum components. The atomic 
state amplitude grating or matter grating created by the microfabricated grating necessarily implies the existence of 
such coherence. However, for L <C Lt, this coherence is not evidence for quantum scattering since it corresponds 
simply to the classical shadow effect. The relative phase associated with the coherence associated with two transverse 
momentum components, p and p + Hq, with hq = 2irh/d, can be written as 

4> = q — t + UJ„t, 

m 

where t = L/uq and tu q = hq 2 /2m is a recoil frequency. The first term is a Doppler shift, independent h - it has a 
classical origin. The second term is a recoil shift which vanishes in the limit h ~ 0. It is a purely quantum term. 
Thus, for Lj q t <C 1, which is equivalent to the condition L <C Lt, the scattering is classical and the coherence between 
transverse momentum components corresponds to a classical Doppler shift. It should now be clear that the fringes 
produced in the two-grating interferometer have nothing to do with quantization of the atomic motion when L <C Lt- 
They result simply from a shadow effect and constitute a moire pattern. For grating separations L > Lt, quantum 
diffraction effects begin to play a role; however, the conditions for fringe production are still determined by classical 
considerations related to cancellation of a Doppler phase. 

We return now to atom interferometers which use optical fields to create atomic coherence between internal states, 
external states, or both. Depending on which coherence is created and probed, one can classify such atom interferom- 
eters into three general categories. Classical atom interferometers (CAI) in no way rely on quantization of the atoms' 
center-of-mass motion for their operation. The CAI operate in the regimes u q t <C 1 or L < Lt, for which quantum 
scattering effects can be ignored. If one restricts the use of the term "atom interferometer" to those devices which 
depend on quantization of the atoms' center-of-mass motion for their operation, the CAI would not qualify. Examples 
of CAI include those based on radio-frequency |ll]] or optical [|l2|-|l4[ Ramsey fringes. 

Matter wave atom interferometers operate only in the regime uj q t > I or L > Lt- We distinguish matter wave 
atom interferometers in which the total atomic density is monitored from those in which an internal state coherence 
or population is probed. The former we refer to simply as matter wave atom interferometers (MWAI) and the latter 
as internal state matter wave atom interferometers (ISMWAI). An important feature which distinguishes MWAI from 
ISMWAI is the dependence of the observed signal on the parameter ui q t or Lj Lt- The ISMWAI signal does not vanish 
in the limit uj q t C 1 or L < Lt, but does vanish in this limit for MWAI. In other words, MWAI depend critically 
on the quantization of the atoms' center of mass motion for their operation, whereas ISMWAI can operate in the 
classical or quantum scattering limits. Although matter wave effects can modify the signals in ISMWAI, many features 
of the signals are determined by classical considerations, as in CAI. Examples of ISMWAI interferometers include 
those which explore recoil splitting of optical Ramsey fringes in the frequency domain p^5Hlq,p 191, and time-domain 
"121 



Matter wave atom interferometers can be constructed using off- resonant standing wave optical fields. Such fields 
act as phase gratings for the matter waves. Atoms enter a field interaction region in their ground state and leave 
the region with their spatial density unaffected. On the other hand, the phase of the ground state wave function is 
spatially modulated, corresponding to a coherence between momentum states differing by nhq (n is an integer and 
4ir/q is the period of the standing wave field). As the matter wave evolves freely following the field interaction region, 
the phase modulated wave function is transformed into a spatially modulated density, which can be probed using 
various techniques. Atom interferometers can be constructed using two or more field interaction zones |22|-|24|]. In 
the transient echo experiment of Cahn et al. |23f| , the vanishing of spatial modulation at the echo times t = 2T and 
t = 3T (but not in the immediate vicinity of the echo times) is consistent with MWAI theory | p5[ . Cahn et al. also 
found that the spatial modulation vanished in the classical scattering region, uj q T -C 1, as it must for any MWAI. 

It is also possible to create MWAI by resonantly driving one-photon transitions between electronic states [^5| or 
two-photon Raman transitions between ground state sublevels using optical fields. The lifetime of the electronic or 
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ground state coherence must be greater than w" 1 , a condition that restricts one to intercombination lines in the case 
of electronic state coherence. It is possible to show that the sum of the spatially modulated atomic gratings 
created by the fields in each of the states involved in the transition (that is, the total atomic density) vanishes in the 
classical scattering limit, but is nonvanishing in the quantum scattering domain, uo q t > 1. This effect has not yet been 
observed. If one were to probe the total atomic density rather than some internal state density in the experiments of 
|^6| , p7|Jl9| , the sought after effect should be observable. 

In this paper, we describe a new type of MWAI [^6| that also uses resonant rather than off-resonant fields. Two-level 
atoms pass through a standing wave, optical field that resonantly drives a closed two-level transition in the atoms. 
In contrast to the scheme described above, it is now assumed that the excited state lifetime is much less than oj" 1 . 
As a result of the atom-field interaction, an excited state matter grating in the excited state is created, as is a "hole 
grating" in the ground state. If one neglects all effects related to quantization of the atoms' center-of-mass motion, 
it is easy to show that, after the excited state decay back to the ground state, the ground state density is uniform. 
In some sense, the excited state grating fills the "hole" in the ground state density that had been created by the 
field. When the recoil associated with stimulated ||^] and spontaneous [^8 29 processes are included, however, new 
features appear in the ground state density matrix. It is not possible to interpret the interaction in simple terms 
as was the case for a phase grating, but the net result is that recoil leads to a nonvanishing spatially modulated 
atomic density for uj q t > 1 or L > Lt- In analyzing the signal, it will prove useful to separate the contributions from 
spontaneous and stimulated processes. Immediately following the excited state decay, these two contributions cancel 
one another, but as time progresses, the contribution from spontaneous processes disappears, leaving a net modulated 
ground state density. When two field interaction zones are used, echo-like phenomena can occur. It will be seen that 
both stimulated and spontaneous processes contribute to the echo signals, even for pulse separations T ^S> wJ 1 . 

The paper is arranged as follows. The change in the atomic density matrix following the interaction of atoms with a 
standing wave, optical field is calculated in Sec. II. In Sec. Ill we consider atom interference using a single atom-field 
interaction zone and a homogeneous distribution of atomic velocities. It is seen that focussing of the atoms, similar 
to that found with phase gratings, can also occur using resonant fields. The conditions under which the spatially 
modulated atomic density becomes a periodic function of t (Talbot effect) arc established. The atomic density following 
two atom-field interaction zones is calculated for an inhomogeneously broadened transition (Talbot-Lau regime) in 
Sec. IV. The results are summarized in Sec. V. 



II. BASIC EQUATIONS 

Atom interferometers can operate in the spatial or time domain. In the spatial domain an atomic beam traverses 
one or more field regions. In the time domain, a vapor of cold atoms (or condensate) is subjected to one or more 
radiation pulses. The spatial domain interferometer can be analyzed in the time domain if calculations are carried 
out in the atomic rest frame. Consequently, we restrict our calculations to the time domain. 

Two-level (upper state |e) and ground state \g}) atoms are subjected to two radiation pulses separated by a time 
interval T. Each radiation pulse consists of two traveling wave components having propagation vectors ki and k2, 
respectively, |ki| = |ka| = fe = fi/c, and is the frequency of each field. The total electric field can be written as 

E (r, t) = ee- mt+iC *- r cos (q • r/2) [E m (t) + E 2 g 2 (t - T)] + c.c, (I) 

where Q = (ki + k 2 ) /2, q = (ki — k 2 ) , e is a polarization vector, Ej is the amplitude of pulse j (j — I, 2), and gj (t) 
is a smooth pulse envelope function having width t, centered at t = Tj, with T\ — and T 2 = T. We assume that the 
pulse separation T, pulse duration r, atom- field detuning A = — u> (lu is the atomic transition frequency), upper 
state decay rate T, recoil frequency 

uj q = hq 2 /2m, (2) 

and Doppler shift k,- • v (v is an atomic velocity) satisfy the inequalities 

TT > I, (3a) 

Tt < I, (3b) 

At < I, (3c) 

ki • v < r, (3d) 

u q T > I. (3e) 

Inequality (|3a|) implies that any excited state population created by the fields decays to the ground state in a time 
that is short compared with the time scale of the experiment. Inequalities (|t], Eteh imply that spontaneous decay 
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and atom- field detuning can be neglected during the radiation pulses, while inequality (3d) (atoms cooled below the 
Doppler limit of laser cooling) guarantees that there is negligible Doppler dephasing for times of order of the excited 
state lifetime. Finally condition ( |3e| ) states that we are in the quantum scattering limit. 

Conditions (||) allow one to map out the time development of the density matrix resulting from each radiation pulse 
in three stages: 

1. an impulsive change in the density matrix produced by the atom- field interaction 

2. spontaneous decay of the excited state 

3. free evolution of the density matrix. 

Depending on the specific application, we calculate the atomic density matrix following the first or second pulse. 
Owing to inequalities (||), we can define times before (T^) and following (T^) pulse j such that changes in the 
atomic density resulting from spontaneous decay and atom-field detunings can be neglected in the time interval 
r (+)_ T (-) 

3 3 ' 

A. Stage 1 

During pulse j, the density matrix, in an interaction representation, evolves according to 

ihp=[V,p], (4) 



V = 2TiXj9j (t — Tj) cos (q • r/2) [cos (Q • r) a x — sin (Q • r) <j y ] , (5) 

where \j = —pEj/2h is a Rabi frequency, p is a dipole moment operator matrix element, a x and a v are Pauli matrices, 
and 

1 \ f 





|S> = (I • (6) 



It has been assumed that, during each pulse, any effects arising from atomic motion can be neglected [inequalities ( pb| , 
3d)], which is the reason the kinetic energy term has been omitted from the Hamiltonian in Eq. (0) (Raman-Nath 
approximation). Before the first pulse acts, it is assumed that all atoms are in their ground state. Owing to inequality 
(|3a|), all population is returned to the ground state before the action of the second pulse at t = T 2 . Thus we need only 
calculate the change in the ground and excited state density matrices produced by pulse j, starting from a density 
matrix 



(r,r'; 





Pgg ( r > r '; T ; 



T t ) ) = \ n Lw.rri-A n I. ( 7 ) 




in which all atoms are in their ground state. One can integrate Eq. (|4j) to obtain the density matrix immediately 
following the pulse at time Tj , 

p (r, r'; T, (+) ) = Vj (r) p (r, r'; T^) „j (r') , (8a) 
T]j (r) = cos [ V2#j cos (q • r/2)] — ism [ 1 / 2 9j cos (q • r/2)] [cos (Q • r) a x — sin (Q ■ r) a y ] , (8b) 

where 



A Xj / dt 9j (t) (9) 



is a pulse area. 

It is convenient to use the Wigner representation for the density matrix, 



f dr 

p(r,p,t)= / - —3 exp (-zp ■ r/h) p(r + r/2, r - r/2, t) . 

J (2-kTi) 
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and expand the populations as 

p nn (r, p, t) = ^ pnn (s, P, t) exp [isq • r] . (10) 

3 

Using Eqs. (|8|)-(|Io|) and expanding the sin [ x cos (q • r/2)] and cos [ 1 jiQj cos (q • r/2)] functions in terms of Bessel 
functions, one obtains for the Fourier coefficients 



Pgg 



; < P> T i +) ) = Yj J ^ (^/2) J 2 (e-s') {Oj/2) Pgg [s - s', P ~hq {I - s'/2) , T$ } ] , (11a) 

t,s' 

Pee (s, p, = £ (-l) s ' J 2 *+i (%/2) J 2( ,_ s , )+1 (0,/2) p 99 {s - a', p-ft [q (* - (s' - 1) /2) + Q] , Tj- } } , (lib) 

e,s' 

where J s (x) is a Bessel function of order s. 

B. Stages 2 and 3 

In the next stages of the calculation, the excited state decays to the ground state and the ground state density 
matrix evolves freely following the decay. The Wigner function associated with the excited state population obeys 
the equation of motion 

d \ 

v ■ V Pee (r, p, t) = -T Pee (r, p, t) (12) 



dt 

It then follows from Eq. ( |l2| ) that the s-order Fourier component evolves as 

p ee (s, p, t) = exp [- (r + isq ■ p/m) (t - Tj)] p e e (s, p, (13) 
The excited state repopulates the ground state. The ground state Wigner function is governed by the equation 



^ + v • VJ p gg (r, p, t) = T J dn r N (n r ) p ee (r, p+?ik r , t) , (14) 
from which it follows that the ground state Fourier coefficients evolve as 

+ isq • p/ni^j pgg (s, p,t) = T J dn r N (n r ) p ee (s, p+hk r ,t) (15) 



In these equations, k r is a spontaneous photon wave vector, n r = k r /fc r , and N (n r ) is the normalized probability 
density for the radiation of a photon in the direction n r . 

The solution of Eq. ([l4]) involves both a homogeneous and particular solution which we write as 

Pgg (S, P , t) = P (f) ( S , P , t) + (S, P , t) . (16) 

The homogeneous part, 

pV? (s, p, t) = exp [-isq • p (t - T 3 ) /m] Pgg (s, p, rj +) ) , (17) 

represents the evolution of the ground state in the absence of decay. The S superscript indicates that this part is 
associated purely with stimulated processes. The particular solution, p g ^ (s, p, t) , represents the contribution to the 
ground state density matrix resulting from excitation by the pulse and subsequent decay of the excited state (hence, 
the superscript D for decay). One finds that at times t — Tj ^> the particular solution is given by 

p { ° ] (s,p,t) =exp[-wq-p(t-r i )/m] J dn r N (n r ) (1 + isLj d n q /T)~ l p ee (a, p+hk r , Tj +) ) , (18) 

where 
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LUd = hqk r /m (19) 

and n q = n r -q/g. The quantity su>d is a recoil frequency associated with the spontaneous decay of the sth excited 
state Fourier component. Alternatively, sujd can be viewed as a Doppler shift of the spontaneously emitted photon 
that is dependent on the momentum kick shq which the atom acquires in the excitation process. 

Piecing together Eqs. Jig), (|l7|), (|TT|), one finds that the ground state density matrix for t — Tj ^> is given by 



p gg (s, p, t) = ex P H«l ' P (* ~ T j) M { J 2£ (Oj/2) J 2(e _ al) (9j/2) Pgg [s - a', p-hq (£ - s' /2) , Tj } 



+J21+1 (Oj/2) J 2 (e-s>)+i (Oj/2) J dn r N (n r ) (1 + isuj d n q /T) 

x Pgg [s - s', p-ft (q (I - s'/2) + Q - k r ) , TH] } (20) 

This is the building-block solution which can be used to analyze several possible experimental schemes. 

For u>~^ t — Tj ^s> r^ 1 , recoil effects can be neglected and the terms involving recoil momenta in the 
arguments of the density matrix elements can be dropped. In that limit the sum over £ of the Bessel functions gives 
S s 'fi, and Eq. reduces to p gg (s, p, t) = cxp [— isq • p (t — Tj) /m] p gg (s, p, T^), which implies that 

Pgg (r, P ,t) - Pgg (t- ^ (t - Tj ) , P^H ) . (21) 

As expected, the density matrix simply undergoes a classical translation if recoil is neglected. 



III. ONE FIELD-INTERACTION ZONE 



When an electromagnetic wave passes through an amplitude or phase grating, the diffraction pattern as a function 
of the distance from the grating is a periodic function of the Talbot length, or, equivalently, of uj q t. Similar effects 
occur for matter waves and have been observed experimentally by Chapman et al |30| and Nowak et al |3l| . In these 
experiments, one sent ground state atoms through a microfabricated grating. When atoms are excited by resonant, 
optical fields, the diffraction pattern is no longer strictly periodic owing to spontaneous decay. In this section we 
calculate the atomic density following the interaction of a highly collimated atomic beam (quO^t <C 1) or a condensate 
(qut <C 1) with a resonant optical field. In the atomic rest frame, the field appears as a pulse centered at Ti = 0. The 

initial Wigner distribution is taken as p gg ^r, p,T^ ^ = W (p) p2] j, corresponding to Fourier components 

p gg (s^T^) = 6 S , W (p) . (22) 

In this equation, W (p) is the momentum distribution of the atoms. 

It is convenient to orient the x, y, z axes along the mutually orthogonal vectors q, Q, and ki x k 2 . It then follows 
from Eqs. (p0[), along with the summation identities 

S^ e) (a) = ^ J2£ (a) Ji{i-v) (a) exp (-Ma) 

e 

= i cxp~ WQ / 2 { J 2u [2a cos (a/4)] + (-1)" J 2l/ [2a sin (a/4)]} , (23a) 

S i° ] («) = J! J2l+1 J 2(<^) + 1 («) eX P 
£ 

= -\ exp-'^- 1 ^ 2 { J 2v [2a cos (a/4)] - (-1)" J 2 u [2a sin (a/4)]} , (23b) 
that the Fourier component at time t, 

Pgg (S, t)= dppgg (s, P,t) , (24) 

is given by 
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Pgg ( s > *) = Va ( ex P H«q ' P*A»)) { Jis [Oi sin (</> Ts (st) /2)] [1 + C (0t d (at) , «V r )] 

+(-l) s J 2s [0i cos (0 Ts (ai) /2)] [1 - C {^To (st) , scj d /T)}} , (25) 

where (. . .) represents an average over atomic momenta, 

<j> Ta (*) = UJqt (26) 

is the Talbot phase associated with stimulated processes, 

<t>T D (*) = w d t (27) 
is the Talbot phase associated with spontaneous processes, and 

C (a, (i)= [ dnN (n) exp (ian x ) (1 + ifln x y l . (28) 



Terms involving C (a, (3) in Eq. ( |25| ) are connected with spontaneous processes while the remaining terms arise from 
stimulated processes. 

It follows from Eq. (p5|) that the spatially homogeneous part of the atomic density is unchanged, 

Pgg(0,t) = l, (29) 

consistent with probability conservation in the closed two-level system. We are interested primarily in the time 
dependence of the Fourier components having s ^ 0, since these components determine the spatial modulation of the 
atomic density. The maximum value of s entering the summation in Eq. ( ^fj| ) is of order max. {1,8%}. It is assumed 
in this section that the Doppler broadening is small, 

max{l,0i}q- pt/m < 1, (30) 



which allows one to set (exp (— isq • pt/m)} equal to unity in Eq. ([25|). Note that, even without this factor, the general 
expression (^5]) is not a periodic function of time. Owing to spontaneous emission the Talbot effect is destroyed. We 
will see below that, for sufficiently large times, the Talbot periodicity is restored. 

If both traveling wave components of the field are linearly-polarized along z and if the ground state angular 
momentum is equal to 0, then 

JV(n)= 3 /s„ (1-nl) . (31) 
We assume that the recoil frequency is sufficiently small to ensure that 

max{l,0i}<Vr < 1. (32) 
As a consequence one need only evaluate (|8|) at = 0. Using the identity Q 

J dnriiUk exp[ia • n] = 47ra~ 3 |(sina — a cosa) S^h H — -j— [(a 2 — 3) sina + 3acosa] | (33) 

one finds 

C (a, 0) = 3 / 2 a -3 [a cos a + (a 2 - l) sin a] . (34) 

Note that C (a, 0) — 1 - 3a 2 /16 for a < 1 and C (a, 0) — ^ sina for a > 1. 

The time scale of the transient response (E5j) is determined by the recoil frequencies, 



Lu d = 4w fe sin (6/2) , u> q — Auj k sin 2 (9/2) (35a) 

where is the angle between the wave vectors and k 2 . There are essentially two time scales in the problem, one 
associated with spontaneous processes, Td — uj^ 1 , and one associated with stimulated processes T q = uj" 1 . For Fourier 
components having s =^ 0, it is possible to isolate the contribution from stimulated processes since the contribution 
from spontaneous processes becomes negligible for t 3> r d . It is not difficult to understand why spontaneous processes 
contribute a negligible amount in this limit. The recoil phase factor associated with spontaneous processes |$4j is 
e ik r -qt/2m _ e in rx u d t i As mentioned above, this phase factor can be viewed as a recoil related Doppler phase. When 
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summed over all directions of the spontaneously emitted photon, it averages to zero for u>dt ^> 1. As a consequence, 
one finds from Eq. (^BJ) that, for ojdt 3> 1 and s ^ 0, 

p gg (*, t) = l / 2 { J 2s [0! sin (0 Ts (st) /2)] + (-1) s J 2s fa cos (0 Ts {st) /2)]} . (36) 



In effect, Eq. (|36|) represents a periodic rephasing (Talbot effect) of the ground state amplitude grating that was 
created immediately following the radiation pulse, since for u>dt 3> 1 the spontaneous contribution no longer plays a 
role. One could have equally well ionized all the excited state atoms immediately following the pulse. The Fourier 
components, p gg (s,t) , and total density, p gg (r, t) = Y"l p„ a (s,t) e lsq r , are periodic functions having period equal to 
lnjuiq and can be used to measure recoil frequency |23| , |35| , but, in contrast to scattering by phase gratings, there is 
no time for which the density is uniform when ujt 3> 1. When all Fourier components are taken into account, and for 
large pulse areas (but not so large as to violate the Raman-Nath approximation), 

(w„t) -1 » e x » 1, (37) 

it can be shown that the atoms are focused by the field. This new regime of atom focusing, as well as its relation to 
focusing by phase gratings, will be considered in a future publication. 

For earlier times, when t < uj^ 1 < LOq 1 ^ the spontaneous term contributes to the atomic density. If t <C w^ 1 , the 
total density is approximately uniform since spontaneous decay " refills" the " hole" in the ground state that is created 
by the radiation pulse. From Eqs. (|25|, HJ), one finds that 

P 99 (*, t) * ^ (s"atf + (~1) S i^dtf J 2s (6i) ■ (38) 
The Fourier components build up as t 2 when t <C w^ 1 - When the angle between wave vectors is small (O < 1) , and 



for somewhat larger times, t ~ uJ d <C uj q 1 , Eqs. (25) reduces to 



p gg (s,t)= 1 /2(-l)V 2s [0 1 ][l-C(0 TD (st),O)], (39) 

allowing one to isolate the contribution from spontaneous processes. The Fourier component p gg (l,i) is plotted in 
Fig. |^ as a function of uidt\ it is not a periodic function of ujdt. The time dependence of Eq. ( |39| ) can serve as a probe 
of the spontaneous emission distribution function [see Eq. (|2^)] . 




For counterpropagating waves (0 = tt), the recoil frequencies uo q and tod coincide and achieve their maximum 
value ojd — uj q — 4wfe. In this limit, atom interference effects occur on the shortest possible time scale. The Fourier 
component p gg (1, t) is plotted as a function of uj q t in Fig. || for 9 = n. Other Fourier components could be shown as 
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well, but p gg (l,t) is the component most easily monitored using backscattering techniques. One sees that, initially, 
p gg (l,t) is aperiodic, but it asymptotically approaches periodic behavior for Udt 3> 1. The area 9\ = 38.9 is that 
for which \p gg (l,t)| achieves its maximum. The area 9i — 4.81 is chosen to maximize the relative contribution of 
spontaneous processes. It was obtained by maximizing the ratio p m /p^, where p m is the maximum of the exact 
expression (p5|), which occurs at Lu q t ~ 1, and is the maximum of the asymptotic expression (|36|) occurring at 

UJqt > 1. 



e,=38.9 




co q t 

FIG. 2. Fourier component p ga (l,t) for a single field interaction zone. The pulse area 6\ — 38.9 (dashed curve) is chosen 
to maximize p gg (l,t), while the area 9i = 4.81 (solid curve) is chosen to maximize the relative contribution from spontaneous 
processes. 



IV. TWO FIELD INTERACTION ZONES 



In the previous section, it was assumed that any Doppler dephasing was negligible on the time scales under con- 
sideration. In this section, we look at the limit in which 

qu^>uj d ,uj q , (40) 

where u is the width of the velocity distribution in the direction of q. To be specific, we take the momentum 
distribution as W(p) = Wx(p±_)W q (p q ), 

W q ( Pq ) = -^—e-PV^) 2 (41) 

is the distribution of momenta in the direction of q and Wj_(Pj_) is the distribution of momenta transverse to q. By 
combining Eqs. ^2^) and (f4l| ) one finds that, following a single pulse, the Fourier components of the density are given 
by 

Pgg (s, t) = e -N«*) 2 /4 J 2s [ 0l(j)Ts ( st ) / 2 ] , (42) 

correct to order uj^/qu. This term survives only for times of order (sgw) -1 , which implies that the argument of the 
Bessel function is of order 9iL0 q /qu. This, in turn, implies that only Fourier components having s ~ max {1, 9iuj q / qu\ 
contribute significantly for times in which Fourier components other than s = are nonvanishing. 

The picture is rather simple. The radiation pulse excites the atoms which then decay back to the ground state 
giving a uniform density. For times t < (squ) , the sth Fourier component begins to build up significantly provided 
9\w q /qu > s, and a modulated atomic density appears. For times t > (qu) , all spatial modulation has been 
washed out as a result of Doppler dephasing. Since this time is shorter than the inverse recoil frequencies, the main 
features of the time dependence p gg (s,t) found in Sec. Ill for homogeneous broadening never can be realized in this 
inhomogeneously broadened sample. 
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The spatial modulation is not lost, however, and can be restored using echo techniques if a second pulse is applied 
at some time T following the initial pulse. The time dependence of p gg (s, t, T), considered as a function of T, displays 
the same features found in the previous section for p gg (s, t) as a function of t. Specifically, it can be used as a probe 
of spontaneous processes. 

In this section, we consider the atomic response following two pulses centered at t = and t = T 
(Ti = and T 2 =T). One key point to remember is that the relevant time window for which the modulation is 
restored is of order (squ) . Thus, it is possible that the Doppler dephasing associated with spontaneous emission 
no longer plays a critical role in killing off the Fourier components, since this dephasing is negligible on a time 
scale (squ) . We shall see this to be the case; as a consequence the Fourier components at the "echo times" have 
contributions from both the stimulated and spontaneous terms, even for uj^t 3> 1. 

The Fourier components of the density at times t — T 2 S> following the second pulse can be written as 

Pgg (a, p,t) = pf g S) (a, p,t) + Pg D g S) {a, p,t) + P f g D) (s, p,i) + Pg D g D] (a, p,t) , (43) 

where pgg (s,p,t) represents the contribution from stimulated [K = S) or spontaneous (K = D) processes following 
the second pulse that depend on the stimulated (/ = S) or spontaneous (I — D) component of the ground state 
densit y m atrix Fourier components that were created by the first pulse. Consider, for example p g S g S ^ (a, p,t) . Using 
Eqs. (11a, [I?], ^l|), one finds that the Fourier component 

P { 9 S 9 S) (s,t) = J dpp(f>(*,p,t), (44) 

is given by 

4f ) (*< *) = E exp {-« w - T ) + ( s - s ') T f / 4 } 

s' 

x exp{-2io; 9 [h -(a- a') /2] [a(t -T) + (a- s')T}} J Ml (8J2) J 2{h - S+S ') (0i/2) 
ii 

x ex P [-^co q (t - T) (£ 2 - a' /2)\ J 2l2 (0 2 /2) J 2( ^_ s0 (0 2 /2) (45) 

When the time separation between pulses is larger than the inverse Doppler width, 

guT»l, (46) 

the average over momenta leads to a non-vanishing contribution only if 

a(t-T) + {a-a')T<(qu)-\ (47) 

Inequality ( [47| ) can be satisfied in the vicinity of the echo times, t e , defined as 

n' , s 

t e = -T, 48 
n 



where n! and n (n! > n) are positive integers having no common factors, provided that 



(49) 



where r is an integer. Setting 

n' 

St = t - t e = t T, 

n 

and using inequality (^7|) and Eqs. (p9j) . one finds that the Doppler phase is nondestructive for times 

St < 1/ (nrqu) 

Although not indicated explicitly, St is a function of t, n, n' and T. 

Since nr = s, Eq. ( |49"| ) implies that the (nr)th Fourier component (r — 1, 2, 3 . . .) is nonvanishing in the vicinity of 
the echo time. For example, if n = 1, all the Fourier components contribute near the echo times t = n'T (n' = 2, 3 . . .), 
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corresponding to a macroscopic atomic grating having period A/ [2 sin (0/2)]; if n = 2, the (2r)th Fourier components 
contribute near the echo times t = n'T/2 (n! = 3,5,7. . .), corresponding to a macroscopic atomic grating having 
period A/ [4 sin (0/2)1. In this manner one can generate macroscopic atomic gratings having period A/ [2n sin (0/2)]. 
Note that condition (M) for a nondestructive Doppler phase is a classical condition since it does not contain h. The 
shape of the echo signal about the echo times and the dependence of the signal on the time separation of the pulses is 
determined by effects related to quantum scattering, but the actual location of the signals is determined by classical 
consideration s onl y [||. 

Using Eq. (23a) for the sums over l\ and i 2: one finds that, at St = t — (n'/n)T, 

pff* (nr,5t) = 1 / 4 exp[-(7irgu^) 2 /4] 

X { J 2r (n'-n) [9l Sin {<hs ( nr5t ) / 2 )1 + (- 1 ) 1 ™ J2r(n'-n) [#1 C «S (<f>T s (nrSt) /2)]| 

x { J 2rn > [02 sin (r0TL s /2)] + (-l) nV J 2rn > [O2 cos (r0 TL /2)]} , (50) 
where 4>tl s is a Talbot-Lau phase, defined as 

4>tl s = (n - n)uj q T. (51) 

It is easy to show that the Fourier component p g ^ (nr, St — 0) does not vanish for s 7^ 0. We shall return to this 
point shortly. 

The remaining terms in Eq. (43) can be evaluated in the same manner, and one finds 

p gg (nr,St,T) =^(St)F(T), (52a) 
$ (St) = V2 exp [-(nrquStf/4\ { J 2r ( n >-n) [0i sin (cb Ts (nrSt) /2)] [1 + C (0 Td (nrSt) , -r (n' - n) cj d /T)\ 

+ (-l) r ( n '- n h 2r{n ,_ n) [0! cos (0 Ts (nrSt) /2)] [1 - C (ct> Tn (nrSt) , -r (n' - n) uj d /T)]} , (52b) 
F(T) - 1 / 2 {J2m> [62Sm(r<j) TLs /2)][l + C (r0 TLD ,rnLu q /T)] 

+ (-1)™ J 2rn ' [^2 cos (rc/) TLs /2)] [1 - C (r4>TL D , rnu) q A 1 )]}, (52c) 

where 



(n' — n) u>dT. (53) 



The result (52a) is the product of a term, <I> (St) , giving the time dependence of the Fourier component in the vicinity 
of the echo time, and a term, F (T) , giving its dependence on the time separation of the pul ses. 

We consider these terms separately, starting with $ (St). The exponential factor in Eq. (52b) leads to a nonvanishing 
contribution to the atomic density only for times nrSt ~ (qu)^ 1 <C u^ 1 < ^q 1 ■ In this limit, and for LUd/T <C 1, Eq. 
(52t) reduces to 

$ (St) = exp [~(nrqu8tf/A] J 2r(n ,_ n) [#i0t s (rotft) /2] . (54) 

Comparing this expression with Eq. (^), one sees that the atomic density near the echo times mirrors the atomic 
density in the time interval St following the first pulse. Thus the dependence of the Fourier component s — nr near 
the echo time can be understood in terms of the dependence of the sth Fourier component at a time St following the 
first pulse. In this time interval, only those Fourier components having s < 6iuj q /qu are created with nonnegligible 
amplitude. 

It is important to note that the recoil dephasing responsible for the "washing out" of the spontaneous contribution 
in the Talbot effect plays no role here, since it is negligibly small in the time interval St ~ (qu)^ 1 . The recoil dephasing 
during the spo ntaneous decay of the excited state does provide a small correction to Eq. (|54|) for St = 0. It follows 
from Eq. ( [52b| ) that, to lowest order in LJd/T <C 1 

$ (St = 0) « 7s (-l) r ( n '- n ) [r (n 1 - n) u d /T] 2 J 2r(n ,_ n) (9,) . (55) 

We return now to the dependence of the Fourier components on T, given by F (T) . The echo configuration considered 
in this section is the same as that which leads to the Talbot-Lau effect. In the (matter wave) Talbot-Lau effect, the 
atomic density is a peri odic function of the time separation between the pulses. In contrast to the normal Talbot-Lau 
effect, the density (52a) is not a periodic function of T, owing to the spontaneous contributions to F (T). However, 
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for pulse separations T ^> uj d , the spontaneous processes no longer contribute to F(T) and one finds the periodic 
dependence 



F (T) « Va { J 2 rn< [to sin (r0 Tis /2)] + (-1)™' J 2r „, [0 2 cos (r^ TLs /2)]} , 
which is reminiscent of Eq. (|36"|). For shorter pulse separations, T -C w^" 1 , the function F(T) builds up as 

F (T) « i (-1)™' J 2r „, (0 a ) [r (n' - n) ^T] 2 . 
[compare with Eq. (38)]. When the angle between ki and k 2 is small, one finds that, for T ~ wj 1 <C w^ 1 , 

F (T) = V 2 (-1)" V J Wr [6 2 ] [l-C {t^ TLd , 0)] 



(56) 



(57) 



(58) 



which is to be compared with Eq. (p9). The Talbot-Lau dependence F(T) [See Eq. (52c)] is qualitatively similar to 



the Talbot dependence of p gg (s,t) [See Eq. (|25j)]. In the limit that 8< 1, and for T 
[Eq. (p9|)], the agreement is quantitative. 



Eq. 
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FIG. 3. A graph of the function F (T) that gives the dependence of the Fourier component p gg (l, 5t,T) on the separation 
between pulses. The pulse area 82 = 7.52 is chosen to maximize F(T), while the area #2 = 10.75 is chosen to maximize the 
relative contribution from spontaneous processes. 



For the case of counterpropagating waves (0 = n) , the dependence F (T) is plotted in Fig. || for n = 1, r = 1 
(s = nr = 1) and n' = 2. The pulse area (9 2 = 7.52 was chosen to maximize F(T), while the area 82 = 10.75 was 
chosen to maximize the contribution from spontaneous processes. The value 62 = 10.75 was obtained by maximizing 
the ratio F m /F^, where F m is the maximum of the exact expression ( |52c| ), which occurs at uj q T ~ 1, and F^ is the 
maximum of the asymptotic expression (|56| ) occurring at Lo q T ^> 1. 



V. SUMMARY 



We have described a new type of matter wave atom interferometer (MWAI). One or two standing wave, resonant 
pulses interact with an ensemble of atoms. Atomic motion during the pulse is neglected (Raman-Nath approximation), 
as is spontaneous emission. As a result of the atom- field interactions, the total atomic density acquires a spatial 
modulation that can be attributable solely to matter-wave interference - the signals arise only for times greater than 
or comparable with the inverse recoil frequency uj^ 1 or uj^ 1 . Spontaneous processes destroy the periodicity of the 
Talbot or Talbot-Lau signals. However, for sufficiently long times, the spatially modulated atomic density becomes a 
periodic function of oj q t (Talbot effect) or cu q T (Talbot-Lau effect). 

To observe the Talbot effect discussed in Sec. Ill, one can use a highly collimated (qu8b <§; w q ) atomic beam that 
is sent through a field interaction region. The fields can be pulsed, if necessary, to ensure that the interaction time is 
much less than the excited state lifetime. The modulated atomic density can be monitored by scattering a probe field 
off the atoms or by directly depositing the atoms on a substrate. One might also contemplate doing this experiment 
in the time domain, using a Bose condensate. The Talbot-Lau effect can be observed either in the spatial domain 
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(using an atomic beam having an appropriate angular divergence) or in the time domain, using a laser cooled and 
trapped vapor. 

Finally, we would like to comment on the fact that the Talbot-Lau Fourier components do not vanish identically 
for St — [see Eq. (p^)]. The amplitude of these components is of order (uid/T) 2 , reflecting the contribution of recoil 
dephasing on the time scale of the excited state lifetime. Although we did not give the equation in the text, there is 
also a contribution to the Talbot Fourier components of order (u)d/T) 2 . These contributions reflect the "opening" of 
the closed two-level system by the recoil associated with spontaneous emission. As in the recoil-induced resonances 
p6| , the opening of the system is connected with quantum scattering - it vanishes in the limit that % ~ 0. 

This situation differs from that involving phase gratings on open, two-level transitions. Imagine that the atoms have 
two ground states, g and g' , to which the excited state can decay, but that the field drives only the g — e transition. By 
using a far detuned field, spontaneous emission to state <?' can be suppressed by a factor (r/A) 2 . Following decay, but 
for times much less than the inverse recoil time, the population density p gg is spatially modulated to order (r/A) 2 , as 
is pg'g', but the total density (p gg + p g ' g ') is uniform. The opening of e — g, two-level system in this case has nothing 
to do with quantum scattering effects. 
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